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Abstract. Analyzing IUE ultraviolet spectra of j3 Cep pulsating stars we have noticed that multiperiodic variables 
have a larger mean metal abundance in the photosphere, [m/H], than monoperiodic ones. We apply statistical 
tests to verify this dichotomy. We obtain that, with a large probability, the multiperiodic (3 Cep stars have greater 
values of [m/H]. This result is consistent with the linear non-adiabatic theory of pulsation of early B-type stars. 
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1. Introduction 

The (3 Cephei pulsators are a well-known and well studied 
group of early-type pulsating stars. Oscillations in these 
stars are strictly connected with the metal abundance Z, 
as they are driven by the classical K-mechanism operating 
in the layer of the metal opacity bump (T « 2 • 10 5 K) 
caused by the lines of Fe-group elements (Moskalik & 
Dziembowski 1992, Dziembowski & Pamyatnykh 1993, 
Kiriakidis et al. 1992 and Gautschy & Saio 1993). The 
size of the instability domain of (3 Cep stars is very sen- 
sitive to the heavy element abundance Z 7 and is smaller 
for lower values of Z (Dziembowski & Pamyatnykh 1993, 
Pamyatnykh 1999). Therefore, the information about the 
metal abundance in (3 Cep stars is of special importance. 

In previous works (Daszyhska 2001, Daszyhska et al. 
2003) we derived the metallicity parameter in the atmo- 
sphere, [m/H], and the mean stellar parameters for all (3 
Cephei stars monitored by the International Ultraviolet 
Explorer. We showed that the metallicity values obtained 
are independent of effective temperature and surface grav- 
ity. Additionally, they are not correlated with the stel- 
lar rotation and any of the pulsational parameters, like a 
dominant period or amplitudes of the light and the radial 
velocity variations. 

In the scope of her PhD thesis, Daszyhska (2001) sug- 
gested that, in general, the multiperiodic (3 Cep stars have 
higher values of [m/H] than the monoperiodic ones. We 
mentioned this result in Daszyhska et al. (2003), and we 
confirmed it in Niemczura & Daszynska-Daszkiewicz 2004 
(Paper I), where we applied a more accurate method of 
determination of stellar parameters. Using the results of 
Paper I, we verify here this hypothesis by means of statis- 
tical tests. 



The paper is composed as follows. In Sect. 2 we remind 
derived values of [m/H]. Sect. 3 contains the test of statis- 
tical hypotheses about equality of two means. In Sect. 4 
we discuss the dichotomy of f3 Cep variables with regard 
to the metallicity. Conclusions are given in Sect. 5. 

2. Metallicity of (3 Cephei variables from UV 
spectra 

In previous work (Paper I) we obtained the values of 
metallicity in the atmosphere for 31 field (3 Cep stars and 
16 ones belonging to three open clusters: NGC 3293, NGC 
4755 and NGC 6231. Our approach consisted of fitting 
theoretical fluxes of Kurucz (1996) to low resolution IUE 
spectra, using the method based on the least-squares op- 
timalization algorithm (Bevington 1969). Errors were es- 
timated by the method of bootstrap resampling. For more 
details about the method we refer the reader to Niemczura 
(2003). 

We analyzed IUE spectra from two calibrations: 
NEWSIPS and INES. The estimated values of the sam- 
ple mean, standard deviation and standard error of the 
mean, for [m/H] derived from these two date sets, are 

NEWSIPS: 



[m/H] = -0.134, S[m/H] = 0.201, = 0.029 



INES: 



[m/H] = -0.125, s [m/H ] = 0.203, s- 



m/H] 



0.030. 



Send offprint requests to: J. Daszynska-Daszkiewicz 



As we have shown in Paper I, these two populations fol- 
low a normal distribution. Applying the two sample t- 
Tcst and the one-way ANOVA test, we found that the 
two means are not significantly different up to the signif- 
icance level a — 0.83. Also, these two data sets give the 
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Table 1. The values of [m/H] for the field Cep stars. 



HD 


Name 


P[d] 


[m/H] 


Multi. 


886 


7 Peg 


0.15175 


-0.04 ±0.08 


N 


16582 


S Cet 


0.16114 


-0.24 ±0.09 


N 


29248 


v Eri 


0.17351 


0.05 ± 0.09 


Y 


44743 


P CMa 


0.25002 


0.04 ±0.10 


Y 


46328 


CMa 


0.20958 


-0.33 ±0.19 


N 


50707 


15 CMa 


0.18456 


0.05 ±0.10 


Y 


52918 


19 Mon 


0.1912 


0.15 ±0.10 


Y 


56014 


27 CMa 


0.0918 


-0.47 ±0.18 


N 


111123 


P Cru 


0.19120 


-0.14 ±0.13 


Y 


116658 


a Vir 


0.1738 


-0.12 ±0.08 


N 


118716 


e Cen 


0.16961 


-0.14 ±0.10 


Y 


120307 


v Cen 


0.195 


-0.26 ±0.10 


Y 


122451 


P Cen 


0.157 


-0.03 ±0.15 


Y 


126341 


r 1 Lup 


0.17735 


-0.39 ±0.16 


N 


129056 


a Lup 


0.25985 


0.04 ±0.10 


Y 


129557 


BU Cir 


0.12755 


-0.27 ±0.11 


Y 


129929 


V836 Cen 


0.14313 


-0.05 ±0.10 


Y 


132058 


P Lup 


0.232 


-0.35 ±0.11 


Y 


136298 


5 Lup 


0.1982 


-0.25 ±0.11 


N 


144470 


LU 1 SCO 


0.067 


-0.45 ±0.15 


N 


147165 


a Sco 


0.24684 


-0.20 ±0.20 


Y 


157056 


9 Oph 


0.14053 


-0.15 ±0.12 


N 


158926 


A Sco 


0.21370 


-0.21 ±0.08 


Y 


160578 


K SCO 


0.19987 


-0.19 ±0.10 


Y 


163472 


V2052 Oph 


0.13989 


-0.25 ±0.16 


N 


199140 


BW Vul 


0.20104 


0.07 ±0.12 


N 


203664 


SY Equ 


0.16587 


-0.01 ±0.21 


N 


205021 


P Cep 


0.19049 


-0.07 ±0.10 


Y 


214993 


12 Lac 


0.19308 


-0.20 ±0.10 


Y 


216916 


16 Lac 


0.16917 


-0.13 ±0.13 


Y 




HN Aqr 


0.15218 


0.21 ±0.09 


N 



same dichotomy of Cep stars with regard to the value of 
metallicity: higher [m/H] for multiperiodic pulsators and 
lower [m/H] for monoperiodic pulsators. Therefore, from 
now on, we will focus only on the NEWSIPS results. 

For the purpose of this work, we remind the obtained 
values of [m/H]. In Table 1 we put the field Cep stars, 
while in Table 2-0 Cep stars belonging to the three open 
clusters. The following columns contain the HD number, 
the star name, the value of the dominant period and the 
value of [m/H]. The information about mono- or multi- 
periodicity at the present knowledge is given in the last 
column (Y=yes, N=no). 

Let us divide the sample of Cep stars in two groups: 
monoperiodic Cep stars (a random variable X) and mul- 
tiperiodic ones (a random variable Y). These samples have 
sizes n = 17 and m = 30, respectively. The sample param- 
eters in these two groups are the following 

x = -0.252, s x = 0.221, s s = 0.054, 

y = -0.068, s y = 0.156, s § = 0.029. 

The same parameters obtained from Gaussian fitting are 

x = -0.258, s x = 0.223, s 2 = 0.056, 



Table 2. The values of [m/H] for the Cep stars in clus- 
ters. 



Name 


ID 


P[d] 


[m/H] 


Multi 


NGC 3293 


V380 Car 


N3293-27 


0.2273 


0.06 ±0.10 


Y 


V400 Car 


N3293-11 


0.1458 


-0.10 ±0.12 


Y 


V401 Car 


N3293-10 


0.1689 


-0.14 ±0.12 


Y 


V404 Car 


N3293-23 


0.1621 


0.27 ±0.07 


Y 


V405 Car 


N3293-14 


0.1524 


0.15 ±0.11 


Y 


V412 Car 


N3293-65 


0.1135 


0.03 ±0.11 


Y 


NGC 4755 


BW Cru 


N4755-F 


0.203 


-0.31 ±0.15 


Y 




N4755-G 


0.156 


-0.56 ±0.20 


N 




N4755-I 


0.179 


-0.29 ±0.10 


N 




N4755-I-13 


0.232 


-0.49 ± 0.12 


N 




N4755-III-01 


0.130 


-0.50 ±0.18 


N 


NGC 6231 


V947 Sco 


N6231-110 


0.1079 


-0.03 ±0.23 


Y 


V920 Sco 


N6231-150 


0.1012 


-0.21 ±0.19 


Y 


V964 Sco 


N6231-238 


0.0878 


0.05 ±0.15 


Y 


V946 Sco 


N6231-261 


0.0988 


0.01 ±0.13 


Y 




N6231-282 


0.1193 


-0.42 ±0.16 


Y 



y = -0.055, s y = 0.188, s v = 0.046. 

Distributions of monoperiodic and multiperiodic stars as 
a function of [m/H], together with Gaussian fits and nor- 
mal curves, are shown in Fig.l. In the left panel, the dis- 
tribution of [m/H] for monoperiodic Cep stars is shown, 
whereas in the right panel the distribution of [m/H] for 
multiperiodic Cep stars is shown. 

Comparing the sample parameters with those derived 
from Gaussian fitting we can say that they are the same 
within the error limits. Additionally we checked the nor- 
mality of these distributions by using the Shapiro- Wilk 
test. Thus, we may assume that X and Y have indepen- 
dent normal distributions N(^x,cr x ) an d N{hy,Oy), re- 
spectively, where /z and a 2 are the mean and the variance, 
respectively. 

3. Statistical tests 

We are interested in comparing the means of the metal- 
licity of the monoperiodic and multiperiodic Cep 
stars, called the random variables X and Y, respectively. 
Because the sample sizes are small, we have to do that in 
two steps. Firstly, we have to test whether the variances of 
these two variables are equal, and then, if we accept this 
hypothesis, we will test the equality of the two means. 

3.1. Test for the variances 

The null hypothesis Ho : o\ = o\ will be tested against 
the alternative hypothesis H\ : o\. For this purpose 
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Fig. 1. The distributions of [m/H] for monoperiodic (3 Cep stars (left panel) and for multiperiodic ones (right panel). 
The normal distribution (solid line) and the Gaussian fit (dotted line) are also shown. 



we use the statistic 



F 



(n-l)S x /(n-l)* x 
(m - 1)5* /(m - 1)4 



C2 

zJL 



which has an F-Snedecore distribution with r\ = n — 1 
and ri = m — 1 degrees of freedom. Thus, F a (n — 1, m — 1) 
is needed to obtain the a significance level. S x and S Y 
are assumed to be unbiased estimates of the correspond- 
ing variances, a\ and a Y , of two independent normal 
distributions with unknown means. To test the hypoth- 
esis Hq : <j x = a Y against a two-sided hypothesis 
Hi : a\ ^ a Y , we assume the critical region given by 



and 



5-2 

> F a /2(n - l,m- 



S 2 

oT ^ F a/2(m - l,n 



1) 



1) 



Taking and as the observed values of unbiased es- 
timates of the variances, we get F = 2.01. At the signif- 
icance level a = 0.05, the critical region, based on the 
random samples of size n — 17 and m = 30, is 



1/4 > ^0.025(16,29) 



and 



si I s z x > F .o2 5 (29, 16) 



or 



Fo.o2 5 (16, 29) < sl/sl < — — = ^0.975(16,29) 

y F . 025(29, 16) 



hence 



2.30 < sl/sl < 0.39 



Because our value of F = 2.01 does not belong to the crit- 
ical region (0,0.39 > U < 2.30, oo), there are no reasons 
to reject the hypothesis H : a\ = a Y - In this case, the 
p- value, i.e. the smallest level of significance, a, at which 
H would be rejected, amounts to 0.83. Thus, with a very 
high probability we may assume the equality of variances. 



3.2. Test for the means 

Given that the hypothesis of equal variances is accepted, 
we are ready to test the hypothesis about the means. If 
random samples of small sizes n and m are taken, the test 
can be based on the statistic 



T = 



X -Y 



Spy/l/n + 1/m 



where 



Sr 



>-l)ff| + (m-l)Sf 
n + m — 2 



The T— statistic has a Student's t distribution with r = 
n+m— 2 degrees of freedom. The critical regions for testing 
the null hypothesis Hq : /ix = Aty against the alternative 
ones are 

(a) \t\>t a / 2 (n + m-2) if H x : [i x ^ [iy , 

(b) t < — t a (n + m — 2) if H 2 : [ix < A*y, 

(c) t > t a {n + m — 2) if H% : fix > fJ-Y- 

We remind that the random variable X represents 
[m/H] for the monoperiodic stars, and the random vari- 
able Y represents [m/H] for the multiperiodic ones. The 
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case (a) means that we test whether the mean values of 
[m/H] for the monoperiodic (3 Cep stars and the multi- 
periodic ones are equal or not. In the case (b) we verify 
the hypothesis of equality of these means against the hy- 
pothesis that the mean for monoperiodic variables is less 
than for the multiperiodic ones. And the last one, (c), is 
the test whether the monoperiodic stars have a greater 
value of [m/H] than the multiperiodic ones. 

If x, y and s^, represent the observed unbiased esti- 
mates of the respective parameters px , Py and o\ = a\ 
of two independent normal distributions, then the empir- 
ical value of the test statistic is t = —3.33. At a — 0.05 
and ?i + to — 2 = 45 we have the following critical regions 
for these three alternative hypotheses 



(a) (-oo,-2.014 > U < 2.014, oo), if H x 

(b) (-oo, -1.680 > if H 2 : p x < Vy, 

(c) < 1.680, oo) if H 3 : p x > Py- 



The values of the observed t belong to the first and the 
second critical regions. Thus Ho is clearly rejected at an 
a = 0.05 significance level in favor of the alternative hy- 
pothesis H i : fix 7^ Py j as well as in favor of the hypothe- 
sis Hi : fix < Py- On the other hand we would not reject 
the hypothesis Hq if the alternative one is H 3 : px > Py ■ 
In this case we have to accept that fix — Py-, but it is 
impossible because H has already been rejected. 

The p— value for the two-side hypothesis, Hi : fix 7^ 
/iy, is 0.002, and for the one-side hypothesis, H2 : px < 
Py, amounts to 0.001. Thus, with a large probability, on 
the basis of these data sets, we can conclude that the mean 
metallicity for multiperiodic (3 Cep stars is higher than 
that for monoperiodic ones. 



4. Discussion of the metallicity dichotomy 

The result from our statistical analysis is in agreement 
with the linear non-adiabatic theory of pulsation of /3 
Cep stars. As was shown by Dziembowski & Pamyatnykh 
(1993), for greater metal abundance we get more unsta- 
ble modes. In Fig. 2 we plot the evolution of oscillation 
frequencies in the Main Sequence phase for four indicated 
masses and two values of the heavy element abundance 
Z. We considered low degree modes with I = 0,1,2. In 
these computations we used the Warsaw-New Jersey stel- 
lar evolution code and the non-adiabatic pulsation code 
of Dziembowski (1977). We took OPAL opacities and did 
not include effects of rotation and convective core over- 
shooting. As we can see from Fig. 2, for the higher Z value 
we get a wider range of unstable frequencies. 

Guided by this property, we plot in Fig. 3 two depen- 
dencies, taking all studied multiperiodic f3 Cep stars. In 
the top panel we show the range of the observed frequen- 
cies as a function of [m/H], whereas in the bottom panel, 
the range normalized by the sum of the lowest and high- 
est frequency as a function of [m/H]. In both cases we 
did not find any correlation. The correlation coefficients, 
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Fig. 2. Frequencies of unstable modes as a function of ef- 
fective temperature in the MS sequence phase of evolution 
for four values of mass: M = 9, 10, 12, 14Mq and two val- 
ues of the heavy element abundance: Z — 0.020 and 0.015. 
Three panels correspond to three values of the spherical 
harmonic degree, I. 



p, calculated including errors as weights, are equal to - 
0.31 and -0.45 for the frequency range and the normalized 
frequency range, respectively. 

This result is not surprising. (3 Cep stars are pul- 
sating variables with small amplitudes both in photom- 
etry as well as in spectroscopy, and we do not know how 
many modes we do observe amongst the excited ones. The 
number of known frequencies increases every year. This 
is mainly thanks to large multi-site and multi-technique 
campaigns, organized already for many [3 Cep variables. 
For example, the last campaign of v Eri (Aerts et al. 2004, 
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Fig. 3. The range of observed frequencies in f3 Cep stars as 
a function of [m/H]. The bottom panel is for the frequency 
range normalized by the sum of the lowest and highest 
frequency. The values of weighted correlation coefficient, 
p, are given in the right-upper corners. 



Handler et al. 2004) gave an increase of the frequency 
range from 0.14 up to 2.28. The data from the 12 & 16 
Lac campaign set up by Gerald Handler are already ana- 
lyzed. Great hopes are pinned also on forthcoming space 
mission, e.g. the COROT project. Therefore the depen- 
dencies in Fig. 3 are just a first trial, which should be 
constantly updated. 



null hypothesis Hq : fix — fJ-Y was rejected in favor of 
Hi : fix 7^ > as well as in favor of H2 : fix < Hy ■ 

The higher value of the mean metal abundance in the 
case of multiperiodic f3 Cep stars is understandable in 
the framework of the linear non-adiabatic theory, which 
predicts more unstable frequencies for higher metal abun- 
dance Z. Therefore, we tried to find a dependence between 
the range of observed frequencies and the metal abun- 
dance. We can not make a statement about any correla- 
tion between these two parameters, because the value of 
the weighted correlation coefficient p amounts to -0.45. In 
fact, we should not expect more, as a number of detected 
modes is much lower than theory predicts. We have to 
remember also that the values of [m/H] give information 
mainly about photospheric metal abundances. Moreover, 
we do not take into account phenomena like diffusion, el- 
ement mixing etc. 

To improve this analysis, a few things must be done. 
Firstly, we need determinations of the [m/H] parameter 
for a larger sample of j3 Cep stars. Secondly, a high ac- 
curacy in the period search analysis is needed to estimate 
better the range of the excited modes in a given star. Here 
put our hope on next multi-site and multi technique cam- 
paigns, as well as on space missions. Finally, the analysis 
of the chemical composition of (3 Cep variables from high 
resolution spectra would be very helpful. 

These results require further studies and must be 
treated with caution. Firstly, they were obtained for small 
samples of stars, and secondly, we do not know how many 
modes from the excited ones have been already identified. 
Our aim was only to show some properties which could 
be extracted from such determinations of the metal abun- 
dance parameter [m/H] for (3 Cep variables. Having in 
mind that from a theoretical point of view all f3 Cep stars 
should be multiperiodic ones, the metallicity dichotomy 
we obtained is rather a qualitative confirmation of the de- 
pendence of the pulsation instability in these stars on the 
heavy element abundance. 

Acknowledgements. We thank Wojtek Dziembowski and 
Alosha Pamyatnykh for the use of the Warsaw-New Jersey evo- 
lutionary code and the non-adiabatic pulsation code. 



5. Conclusions 

Using statistical tests, we have shown that with a very 
high probability the multiperiodic (3 Cephei stars have 
a higher mean metallicity in their atmospheres than the 
monopcriodic ones. The statistical analysis was based on 
our last determinations of the metallicity parameter [m/H] 
(Paper I). Getting a normal distribution of [m/H] for the 
mono- and multiperiodic stars, we applied the test about 
equality of variances of these two samples, represented 
by X and Y random variables. The test based on the 
F statistic allowed clearly to accept the null hypothesis 
Ho : <Ty = u\. Then we could test the hypothesis about 
the mean values of metallicity for the mono- and multi- 
periodic variables. We used the T statistic to define the 
critical regions for a small size of random samples. The 
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